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A  second-order  perturbation  theory  is  developed  for  the  response  of 
slightly  damped  gyroscopic  systems.  The  solution  is  based  on  the  eigen- 
solution  for  undamped  gyroscopic  systems  and  is  expressed  in  terms  of 
real  quantities  alone. 


Introduction 


A  general  theory  exists  for  the  response  of  linear  systems  to 
arbitrary  time-dependent  excitations  (Ref.  1).  The  solution  is  based  on 
the  so-called  transition  matrix.  For  high-order  systems,  the  determination 
of  the  transition  matrix  is  time-consuming,  so  that  the  method  is  not 
particularly  attractive  computationally.  The  situation  is  considerably 
better  for  special  classes  of  systems,  as  shown  in  the  following. 

Undamped  nongyroscopic  systems  are  characterized  by  symmetric  mass 
and  stiffness  matrices.  The  eigensolution  for  such  systems  consists  of 
real  eigenvalues  and  eigenvectors,  where  the  latter  are  orthogonal  with 
respect  to  the  mass  matrix.  Taking  advantage  of  these  properties ,  the 
response  can  be  obtained  without  much  difficulty  (Ref.  1).  If  viscous 
damping  is  present,  then  the  eigensolution  ceases  to  be  real,  even  when 
the  damping  matrix  is  symmetric.  Both  the  eigensolution  and  the  response 


are  significantly  more  difficult  to  obtain  than  for  undamped 
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(Ref.  1).  A  method  for  obtaining  the  eigensolution  and  response  of  un¬ 
damped  gyroscopic  systems  has  been  developed  recently  by  the  first  author 
(Refs.  2  and  3).  Taking  advantage  of  the  fact  that  the  gyroscopic  matrix 
is  skew  symmetric,  the  eigenvalue  problem  can  be  transformed  into  one  in 
terms  of  symmetric  matrices  alone. 

The  problem  of  damped  gyroscopic  systems  is  considerably  more  compli¬ 
cated  than  those  corresponding  to  the  three  special  cases  mentioned  above. 
The  sum  of  the  damping  matrix  and  the  gyroscopic  matrix  is  an  arbitrary 
matrix.  Hence,  any  advantage  resulting  from  the  symmetry  or  skew  symmetry 
of  coefficient  matrices  is  lost,  so  that  one  must  return  to  the  general 
theory. 

This  paper  is  concerned  with  gyroscopic  systems  with  small  damping. 

In  this  case,  damping  can  be  regarded  as  a  perturbation  to  the  undamped 
gyroscopic  system.  Indeed,  a  second-order  perturbation  theory  for  the 
response  of  slightly  damped  gyroscopic  systems  is  developed  in  this  paper. 
The  solution  is  based  on  the  eigensolution  for  undamped  gyroscopic 
systems  and  is  expressed  in  terms  of  real  quantities  alone.  The  theory 
contains  the  case  of  slightly  damped  nongyroscopic  systems  (Ref.  4)  as  a 
special  case. 

2.  Response  of  General  Damped  Gyroscopic  Systems 

The  equations  of  motion  of  a  general  damped  gyroscopic  system  can 
be  written  in  the  matrix  form  (Ref.  1) 

Mq(t)  +  (G  +  C)q(t)  +  Kq(t)  -  Q(t)  (1) 

where  M,  C  and  K  are  n  *  n  real  symmetric  matrices,  G  is  an  n  *  n  real 
skew  symmetric  matrix,  q(t)  is  the  n-dimensional  configuration  vector 
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and  Q(t)  is  the  associated  n-dimensional  force  vector.  We  shall  consider 
the  case  in  which  M  and  K  are  positive  definite. 

The  solution  of  Eq.  (1)  can  be  most  conveniently  derived  by  trans¬ 
forming  the  equation  to  state  form.  To  this  end,  we  introduce  the  2n- 
dimensional  state  vector  and  excitation  vector 

x(t)  -  [qT(t)  j  qT(t)]T 

X(t)  -  CQT(t)  |  OY 

as  well  as  the  2n  *  2n  matrices 

-  ■  [rti] 

k*  - 

where  M*  is  real  symmetric  and  positive  definite  and  K*  is  real  but 
arbitrary.  Definitions  (2)  and  (3)  permit  us  to  rewrite  Eq.  (1)  in  the 
state  form 

M*x(t)  +  K*x(t)  -  X(t)  (4) 

Before  discussing  the  solution  of  Eq.  (4),  we  wish  to  reduce  it  to 
standard  form,  in  which  the  coefficient  matrix  multiplying  x  is  simply  the 
identity  matrix.  Ordinarily,  this  would  mean  premultiplication  of  Eq. 

(4)  by  (M*)  .  Owing  to  the  special  nature  of  the  problem,  however,  there 

exists  a  computationally  superior  procedure.  Indeed,  because  the  matrix 
M*  is  real  symmetric  and  positive  definite,  we  can  use  the  Cholesky 
decomposition  (Ref.  1)  and  write 

M*  -  LLT  (5) 

where  L  is  a  lower  triangular  matrix.  Then,  introducing  the  linear 
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(2a) 

(2b) 

(3a) 

(3b) 


transformation 

LTx(t)  -  u(t)  ,  x(t)  *  L  Tu(t)  (6) 

-T  T  -1  _i  t 

where  L  ■  (L  )  ■  (L  )*,  Eq.  (4)  can  be  reduced  to  the  standard 

form 

u(t)  -  Au(t)  +  U(t)  (7) 

in  which 

A  -  -  L_1K*L-T  (8) 

is  a  real  nonsymmetric  matrix  and 

0(t)  -  L_1X(t)  (9) 

is  a  real  vector. 

The  solution  of  Eq.  (7)  can  be  written  in  the  general  form  involving 
the  convolution  integral  (Ref.  1) 

t 

u(t)  -  4>(t)u(0)  +  /  4>(t,T)U(T)dT  (10) 

0 

where  u(0)  is  the  initial  vector  and 

4(t,r)  -  eA(C'T)  (11) 

is  known  as  the  transition  matrix. 

The  solution  of  Eq.  (7)  can  also  be  obtained  by  modal  analysis, 
which  amounts  to  the  determination  of  the  Jordan  form  for  A.  To  this  end, 
let  us  consider  the  eigenvalue  problem 

Au1  -  Xiyi  ,  i«l,2,...,2n  (12) 

where  and  u^  are  the  eigenvalues  and  eigenvectors  of  A.  For  simplicity, 
we  shall  assume  that  all  the  eigenvalues  are  distinct,  so  that  the  Jordan 
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matrix  Is  diagonal 


A  -  diag[X1]  (13) 

The  eigenvectors  are  known  as  the  right  eigenvectors  cf  A  and  can  be 
arranged  in  the  square  matrix 

U  -  [u1  u2  ...  u2nJ  (14) 


The  adjoint  eigenvalue  problem  is 

-iA  “  Xi-i  »  i“1'2 . 2n  (15) 

and  has  the  same  eigenvalues  X^  as  before,  as  well  as  the  left  eigenvectors 
v^^  (i*l,2,. . . ,2n) .  The  left  eigenvectors  can  be  arranged  in  the  square 
matrix 


[y 


1  ~2 


(16) 


The  right  eigenvectors  are  orthogonal  to  the  left  eigenvectors,  a  property 
known  as  biorthogonality.  The  eigenvectors  can  be  normalized  so  as  to 
satisfy 


T 

VAU  - 


T 

U  V 


21 


(17) 


where  I  is  the  identity  matrix,  in  which  case  the  Jordan  matrix  can  be 
expressed  as 

VTAU  -  2A  (18) 

The  reason  for  the  factor  2  on  the  right  side  of  Eqs.  (17)  and  (18)  will 
become  evident  later. 

Because  the  vectors  u^  on  the  one  hand  and  the  vectors  v^  on  the 
other  hand  are  linearly  independent,  either  set  of  vectors  can  be  taken 
as  a  basis  for  a  linear  vector  space  L^n,  which  implies  that  the  solution 
of  Eq.  (7)  can  be  represented  as  a  linear  combination  of  the  eigenvectors 


i 


r' 


or  of  the  eigenvectors  v^.  Using  the  vectors  as  a  basis,  we  can 
write 

2n 

u(t)  -  Z  u  z  (t)  «  Uz(t)  (19) 

i-1  1 

where  z(t)  is  a  2n-vector  with  components  z^(t).  Introducing  Eq.  (19)  into 

T 

Eq.  (7),  premultiplying  the  result  by  V  and  considering  Eqs.  (17-18),  we 
obtain 

z(t)  -  Az(t)  +  Z(t)  (20) 

where 

Z(t)  -  j  VTU(t)  (21) 

Equation  (20)  represents  a  set  of  2n  independent  equations  having  the 
solutions 

z±(t)  =  eAitZl(0)  +  fQ  eAl(t"T)Z1(T)dT  ,  i«l,2,...,2n  (22) 

or  in  vector  form 

z(t)  -  eAtz(0)  +  /Q  eA(t‘T)Z(T)dx  (23) 

where  we  note  that  eAt  is  a  diagonal  matrix.  The  formal  solution  is 
completed  by  first  introducing  Eq.  (23)  into  Eq.  (19)  and  then  intro¬ 
ducing  the  result  into  the  second  of  Eqs.  (6). 

The  interest  lies  in  systems  with  a  large  number  of  degrees  of 
freedom,  measuring  in  the  hundreds  or  more.  For  such  systems,  the  proce¬ 
dures  described  above  must  receive  closer  scrutiny. 

Solution  (10)  involves  an  integral  containing  the  transition  matrix 
4,  Eq.  (11).  To  generate  the  matrix  4,  it  is  necessary  to  expand  a  power 
series  in  A.  Then,  the  solution  is  obtained  by  performing  the  Indicated 
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integration.  For  large-order  systems,  the  process  can  be  time-consuming. 
In  addition,  computer  roundoff  is  likely  to  introduce  computational 
errors.  The  other  procedure,  namely,  that  based  on  modal  analysis,  is 
also  not  very  attractive  for  high-order  systems.  Indeed,  the  procedure 
requires  the  solution  of  the  algebraic  eigenvalue  problem  for  the  matrix 
A.  But  A,  although  real,  is  not  symmetric.  Hence,  in  general  the  eigen¬ 
values  and  eigenvectors  are  likely  to  be  complex  quantities,  as  opposed 
to  real  quantities  for  real  symmetric  matrices.  Moreover,  computational 
algorithms  for  the  elgensolutlon  of  arbitrary  matrices  are  not  nearly  as 
efficient  as  those  for  real  symmetric  matrices.  Hence,  for  high-order 
matrices,  serious  numerical  difficulties  can  be  encountered,  and  a  dif¬ 
ferent  approach  appears  highly  desirable. 


3.  Undamped  Gyroscopic  Systems 

We  observe  from  Eqs.  (3)  and  (8)  that  in  the  absence  of  damping, 

C  *  0,  the  matrix  A  becomes  skew  symmetric.  In  this  case,  the  eigenvalues 
occur  in  pure  imaginary  complex  conjugate  pairs  and  the  eigenvectors  also 
occur  in  complex  conjugate  pairs,  or 

Au^  *  *  ^i-i  »  i"l»2,...,n  (24a) 

Au^  *  X^^  -  -iu^u^  »  i“l»2,...,n  (24b) 

T 

where  is  the  ^th  natural  frequency  of  oscillation.  Because  A  *  -A, 
the  left  eigenvectors  of  A  are  simply  the  complex  conjugates  u^.  Indeed, 
if  we  write  first 


T  T  T 

VjA  *■  X^v^  ■  iw^v^  *  i"l,2,...,n 

and  then  transpose  the  equations,  we  obtain 

Av^  “  -  X^v^  “  »  i*l,2,...,n 


(25) 


(26) 
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so  that,  comparing  Eqs.  (24b)  and  (26),  we  conclude  that 


v±  *  ui  ,  1*1,2,. . • ,n  (27) 

T 

Hence,  the  2n  eigensolutions  of  A  and  A  can  be  written  in  the  special 
form  of  complex  conjugate  pairings,  or 


X 


i 
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»i 


9 


*i 


»i 
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i* 1,2,. .« ,n 


(28) 


The  above  implies  that  it  is  not  really  necessary  to  solve  the  eigenvalue 
T 

problem  for  A  to  obtain  the  left  eigenvectors  v^. 

Next,  let  us  premultiply  Eqs.  (24)  by  -A  and  write 
2  2 

-A  u^  *  -itUjAu^  *  oi jU  ,  i-l,2,...,n  (29a) 

2—  —  2— 

-A  u^  «  ioi^Au^  -  oi^u^  ,  i*l,2 . n  (29b) 

2  2 
from  which  we  conclude  that  the  matrix  -A  admits  the  same  eigenvalue  oi^ 

—  2 
for  both  eigenvectors  u^  and  u^.  Hence,  the  2n  eigenvalues  of  -A  consist 

2 

of  n  eigenvalues  with  double  multiplicity.  All  the  eigenvalues  of 

2  2 
-A  are  real  and  positive,  which  is  consistent  with  the  fact  that  -A  is 

not  only  real  and  symmetric  but  also  positive  definite.  But  a  real  sym¬ 
metric  matrix  is  known  to  possess  only  real  eigensolutions,  so  that  it  is 
not  necessary  to  work  with  complex  eigenvectors.  Indeed,  because  any 

linear  combination  of  the  eigenvectors  u^  and  u^  is  also  an  eigenvector  of 
2  2 

-A  belonging  to  we  can  simply  choose  the  two  real  eigenvectors 

y  (u^  +  ■  Re  ^  ,  -  y  (u^  -  Uj)  *  Im  -  *  (30) 

In  conclusion,  the  eigenvalue  problem  for  the  undamped  gyroscopic 

system,  Eqs.  (24),  has  been  reduced  to  the  eigenvalue  problem  for  the 

2 

real  symmetric  positive  definite  matrix  -A  .  The  latter  problem  has 
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multiplicity  two,  with  the  real  eigenvectors  ^  and  belonging  to  the 
2 

same  eigenvalue  (i*l,2 , . . . ,n) •  These  are  essentially  the  results  ob¬ 
tained  in  Ref.  2. 

The  reduction  of  the  eigenvalue  problem  for  an  undamped  gyroscopic 

system  to  that  of  a  real  symmetric  positive  definite  matrix  makes  possible 

the  use  of  a  large  variety  of  computationally  efficient  algorithms  for 

the  solution  of  the  latter  problem.  The  question  remains,  however, 

whether  there  is  an  efficient  way  of  producing  the  response  of  a  damped 

gyroscopic  system.  For  arbitrarily  large  damping,  the  answer  must  be 

negative,  as  the  addition  of  the  matrix  C  to  G  destroys  the  skew  symmetry 

2 

of  A  and  the  symmetry  and  positive  definiteness  of  -A  .  On  the  other 
hand,  for  relatively  small  damping,  it  is  possible  to  devise  a  perturba¬ 
tion  scheme  taking  advantage  of  results  obtained  for  the  undamped  case. 
Before  discussing  the  specific  case  of  slightly  damped  gyroscopic 
systems,  we  shall  present  a  second-order  perturbation  theory  for 
arbitrary  real  matrices. 

4.  Second-Order  Perturbation  Theory  for  the  Algebraic  Eigenvalue  Problem 
Let  us  consider  a  2n  x  2n  arbitrary  real  matrix  A.  The  eigenvalue 
problem  for  A  is  given  by  Eq.  (12)  and  the  adjoint  eigenvalue  problem  is 
given  by  Eq.  (15).  The  two  eigenvalue  problems  are  rewritten  as 

Au^  *  »  i*l»2,...,2n  (31a) 

ATvi  •  X1y1  ,  1*1,2 . 2n  (31b) 

where  X^  is  the  ith  eigenvalue  (the  same  for  both  problems)  and  u^  and  v^ 

T 

are  eigenvectors  of  A  and  A  ,  respectively.  As  mentioned  in  Sec.  2,  u^ 
and  v^  are  called  right  and  left  eigenvectors  of  A,  respectively,  and 
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possess  the  biorthogonality  property.  Upon  the  normalization  indicated 
by  Eqs.  (17)  and  (18),  the  biorthonormality  relations  have  the  explicit 


form 


T  T 

*  26ij  >  i,j»l,2,... ,2n 

vTau.  -  u?ATv  *  25  .  ,  i,j-l,2,... ,2n 


(32a) 

(32b) 


i  ~i  ij 

It  should  be  pointed  out  that  eigenvalues  and  eigenvectors  satisfying 
Eqs.  (32)  satisfy  Eqs.  (31)  automatically,  but  the  reverse  is  not  neces¬ 
sarily  true. 

Next,  let  us  assume  that  the  matrix  A  can  be  expressed  as 


A  “  Ao  +  Ax 


(33) 


where  Aq  and  A^  are  known  matrices,  with  the  elements  of  A^  being  of  one 
order  of  magnitude  smaller  than  the  elements  of  A^.  In  fact,  the  matrix 
A  can  be  regarded  as  the  result  of  perturbing  the  matrix  A^  slightly. 
Consequently,  we  shall  refer  to  A^  as  the  unperturbed  matrix,  to  A  as  the 
perturbed  matrix  and  to  A^  as  the  perturbation  matrix.  The  eigenvalue 
problem  associated  with  AQ  will  be  identified  as  the  unperturbed  eigen¬ 
value  problem.  Denoting  the  eigenvalues  of  A^  by  X^,  the  right  eigen¬ 
vectors  by  u^  and  the  left  eigenvectors  by  v^  ,  ana  recognizing  that 
they  all  satisfy  eigenvalue  problems  of  the  type  (31)  and  biorthonormality 
relations  of  the  type  (32) ,  we  can  state  by  analogy  the  unperturbed 
eigenvalue  problem 


Voi  ■ 

X0i*01 

> 

i-1,2,.. . ,2n 

AnU  * 

X  u 

i-1,2,. . . ,2n 

0~0i 

Oi'Oi 

and  the  corresponding  biorthonormality  relations 


(34a) 

(34b) 
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T  T 

~0j~0i  “  -Oj'Oi 


26 


ij 


?0jV0i  "  ^Oi  "  2XOi6ij 


j  • • > 2n 

*  i» j*l»2 , . . . ,2n 


(35a) 

(35b) 


Now  let  us  assume  that  unperturbed  eigenvalues  A  ^  and  eigenvectors 
u^  and  Vq^,  satisfying  Eqs.  (34)  and  (35),  are  known.  Then,  our  objective 
is  to  produce  eigensolutions  of  the  perturbed  eigenvalue  problem  in  terms 
of  the  unperturbed  eigensolutions  and  the  perturbation  matrix  A^.  To  this 
end,  we  seek  solutions  of  the  perturbed  eigenvalue  problem  in  the  form 


xi' 

A0i  +  AU  +  *2i  +  '  ’ ' 

,  1*1,2,. . . ,2n 

(36a) 

»i“ 

~0i  +  ~li  *  ~2i  + 

,  i=l,2,. . . ,2n 

(36b) 

Yi  * 

?0i  +  ?li  +  ?2i  +  •" 

,  1*1,2,. . . ,2n 

(36c) 

where  the  first  subscript  identifies  the  order  of  any  particular  term. 

For  example,  A  is  of  order  one,  and  hence  is  one  order  of  magnitude 

smaller  than  A  and  one  order  of  magnitude  larger  than  A^.  Similarly, 

v^A  is  of  order  two,  and  hence  is  one  order  of  magnitude  smaller  than 
y  and  two  orders  of  magnitude  smaller  than  v^.  We  note  that  the  expan¬ 
sions  (36)  are  not  finite.  Hence,  if  all  terms  in  these  expansions  are 
determined,  these  expansions  must  converge  to  the  eigensolutions  of  the 
perturbed  eigenvalue  problem.  Because  we  shall  determine  and  use  only  a 
finite  number  of  terms,  the  resulting  truncated  expansions  will  provide 
approximate  eigensolutions,  which  must  tend  to  the  unperturbed  eigen¬ 
solutions  as  tends  to  zero.  Furthermore,  as  A^  tends  to  zero  the  terms 
in  any  one  of  Eqs.  (36)  must  maintain  the  same  order  of  magnitude  in 

relation  to  each  other.  For^  example,  we  cannot  have  the  situation  where 

2 

U£^  tends  to  zero  as  A^  and  u^  tends  to  zero  as  A^. 
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The  question  arises  as  to  what  equations  one  should  use  to  determine 


the  perturbations  to  the  eigenvalues  and  eigenvectors.  Lancaster  (Ref. 
4),  Franklin  (Ref.  5)  and  Wilkinson  (Ref.  6)  present  derivations  in  which 
Eqs.  (31)  are  utilized.  Franklin  carries  the  derivation  only  through 
first  order,  while  Lancaster  and  Wilkinson  pursue  the  manipulations 
through  second  order.  Wilkinson  points  out  that  normalization,  as  in 
Eqs.  (32),  is  lost  and  that  the  perturbed  eigenvectors  must  be  renorma¬ 
lized  (Ref.  6,  Sec.  2.7).  To  avoid  this  difficulty,  it  is  suggested  in 
this  paper  that  Eqs.  (32)  be  used  to  determine  the  perturbations  to  the 
unperturbed  eigensolutions.  As  mentioned  previously,  Eqs.  (32)  imply  not 
only  normalization  but  also  satisfaction  of  the  perturbed  eigenvalue 
problem,  Eqs.  (31). 

Substituting  Eqs.  (33)  and  (36)  into  Eqs.  (32)  and  separating 
according  to  order  of  magnitude,  we  obtain 


0(0): 


T 

v  u  *  26 
y0j-0i  ij 


~0jA0'0i  =  2XOi6ij 


j“l*2 , . • • f 2n 


(37a) 


OU):  yjuu  +  yju0i  -  0 

y0jA0'li  +  Y0jAly0i  +  yljA0'0i  "  2Xli6ij 


i,j=l,2, . . . ,2n 


(37b) 


0(2):  ?0j»21  +  ?J»11  +  % *01  ’  ° 

-0jA0~2i  +  ?0jAl*U  +  SljVli  +  » JVoi 


i,j ,=1,2, . . . ,2n  (37c) 


+  y2jV0i  =  2X2i5ij 


As  expected,  Eqs.  (37a)  are  identical  to  Eqs.  (35).  Turning  our 
attention  to  Eqs.  (37b),  we  now  wish  to  determine  u^  and 

Because  u^  is  a  2n-vector  in  the  space  L2n,  and  because  the  vectors 
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Uq^  (i*l,2 , • . . ,2n)  span  this  space,  we  can  represent  as  a  linear  combi¬ 
nation  of  the  Hence  we  write 

2n 

-li",Z®ik-0k  ’  i“1»2»* • • »2n  (38) 


and  we  note  that  the  smallness  of  relative  to  u^  requires  that  the 
coefficients  e  ^  be  small.  Following  the  same  pattern,  we  represent 
v^  as  a  linear  combination  of  the  v  as 


Yli 


2n 
2  Y 
k-1 


ik-Ok 


i*l,2 , . . . ,2n 


(39) 


where  the  coefficients  y„  are  also  small.  Substituting  Eqs.  (38)  and 
(39)  into  Eqs.  (37b)  and  using  Eqs.  (37a),  we  obtain 


6  +  v 

ij  ji 


i ,  j”l , 2, . . . ,n 


X0j€ij  +  X0iYji  *  “  2  '0jAl~0i  Xli5ij 


i,  j*l,2 , . .  •  *2n 


(40a) 

(40b) 


When  i  1  j ,  6^  "  0,  so  that,  solving  for  and  y  we  obtain 


'ij 


-  Y 


y0iAi^0i 


^  2(X0i"X0j) 


i,j-l,2,...,2n,i*j 


(41) 


Alternatively,  for  i  *  j  Eqs.  (40)  become 


*  0  ,  i*l>2, . . . ,2n 

X0i(6ii  +  Yii)  "  “  2  Y0iAl~0i  +  Xli  *  i-1»2,....2n 


(42a) 

(42b) 


Equations  (42)  yield  the  first  order  perturbations  to  the  eigenvalues 


Xli  "  2  Y0iAly0i  ’  i*i»2.---»2n 


(43) 


Moreover,  even  though  Eqs.  (42)  do  not  determine  and  y^  uniquely, 
these  equations  are  clearly  satisfied  if  we  take 
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r 


eii  "  Yii  "  0  ’  1-1,2 . 2n 


Equations  (41),  (43)  and  (44)  determine  the  first-order  corrections 


u^  and  v^  fully.  The  reason  for  our  assumption  that  the  unperturbed 


eigenvalues  be  distinct  is  evident  from  Eqs.  (41).  It  is  easy  to  verify 


from  Eqs.  (41)  and  (43)  that  X^,  u^  and  v^  do  tend  to  zero  as  tends 


to  zero. 


Let  us  now  consider  Eqs.  (37c)  and  obtain  the  second-order  correc¬ 


tions.  Proceeding  as  before,  we  let 


~2i  “  ^ik-Ok  »  1-1,2 . 2n 


(45a) 


?21  ’  ^iTl!r--°k  ’ 


(45b) 


where  and  y^  are  second-order  quantities.  Substituting  Eqs.  (45) 


into  Eqs.  (37c),  and  using  the  biorthonormality  of  the  unperturbed  eigen- 


solutions,  Eqs.  (37a),  we  obtain 


eij  +  J^ik^k  +  Yji  *  °  ’ 


(46a) 


A0jeij  +  2  kf1€iky0jAly0k  +  ^Vik^k  +  2  ^ x Yj ky0kAl~ Oi 


+  \n.y..  -  X- .6 . .  .  .  ,  . 

Oi  ji  2i  ij  ,  i,j-l,2,... ,2n 


(46b) 


Using  Eqs.  (40),  Eqs.  (46)  can  be  reduced  to 


eij  +  Yji  "  ^ Ieikekj  *  i,j-1,2.'**.2n 


(47a) 


X0J€ij  +  X0iYji  =  (Xli  “  Xlj)eij  +  X2i6ij 


Oi  Oj  Ok  lk  kj 


(47b) 


M 


When  i  4  j,  <5^  *  0,  so  that  Eqs.  (47)  can  be  solved  for  6^  and  y^, 
with  the  result 


«  X0i‘X0j 


(Xlj‘Xli)eij  +  kf1(X0k_Xpj)eikekj 


i,j-l,2,...,2n  (48a) 


J1  x0i_x0j  L 


2n 


(Xli_Xlj)6ij  +  kf1(X0i_X0k)eikekj 


i,j-l,2,... ,2n  (48b) 


When  i  *  j,  Eqs.  (47)  become 

'it +  ‘'n  ■  •  1,j'1,2 . 2” 


(49a) 


X0i(6ii  +  Yii) 


X2i  +  J,  (2X0i 

k-1 


X0k)eikeki 


i,j-l,2,... ,2n  (49b) 


Solving  Eqs.  (49)  for  A^,  we  obtain 


2i 


2n 

-  E 

k-1 


(A 


Ok 


-  ‘oi)eik'kl 


i-1,2 . 2n 


(50) 


Although  Eqs.  (49)  do  not  determine  and  y  ^  uniquely,  the  equations 
are  satisfied  if  we  take 


1 

2 


2n 
E  e 
k-1 


ikeki 


i-1,2,. . . ,2n 


(51) 


Equations  (48),  (50)  and  (51)  determine  the  second-order  corrections, 

X2i*  >2i  an<*  ^2i  Recalling  that  A^,  and  y  are  proportional 

2 

to  A^,  we  conclude  that  A^,  u21  and  y2i  are  proportional  to  A^,  as 
anticipated. 
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The  8econd-order  perturbation  solution  derived  in  the  preceding 
section  is  quite  general  and  the  only  restriction  is  that  the  matrix  A 
be  real.  In  the  case  of  slightly  damped  gyroscopic  (or  nongyroscopic) 
systems,  the  matrix  C  can  be  regarded  as  being  of  one  order  of  magnitude 
smaller  than  matrices  G  and  K.  Hence,  recalling  Eqs.  (3b),  (8),  and  (33), 
we  can  write 


(52a) 

(52b) 


One  should  note  that  the  matrix  L  can  be  partitioned,  allowing  a  reduction 

in  computational  effort.  Indeed,  from  Eqs.  (3a)  and  (5),  we  can  write 
L  |  0  " 

L  -  5— +-E-  (53) 

L  i  t-l 

T  T 

where  M  *  and  K  *  L2L2*  Hence»  E<ls*  (52)  can  now  be  written  as 


L  ”1GL ,'T  !  L  "1L 


:SL  -t  _t 

2  Li  : 


-1  -T  i 
L.  aCL.  I  0 

1  1  I 

A,  ■  -  - -t- — 


1  “2 


(54a) 


(54b) 


Note  that  is  a  ^ow,r  triangular  matrix,  so  that  AQ  is  banded 

The  symmetry  and  skew-symmetry  of  and  A^,  respectively,  is  readily 
apparent.  Because  matrix  A^  is  of  one  oro*r  of  magnitude  smaller  than 
Aq,  the  developments  of  the  preceding  section  are  readily  applicable. 

In  fact,  because  matrix  Aq  is  skew  symmetric,  the  unperturbed  eigensolutions 
of  Sec.  4  are  in  reality  the  eigensolutions  of  the  undamped  gyroscopic 


system  discussed  in  Sec.  3.  Recalling  the  results  of  Sec.  4,  we  anticipate 


products  of  the  form  Vq  AjJJq.  Because  only  the  first  n  rows  and  columns 
of  the  2n  x  2n  matrix  are  nonzero,  we  need  use  only  the  upper  halves  of 
vectors  v^  and  u^. 

Let  us  introduce  primed  subscripts  where  a  prime  indicates  that  the 
subscript  in  question  has  been  increased  by  n.  For  example,  j*  =  j  +  n, 
j  *  l,2,...,n.  Using  this  notation  we  have,  from  Eqs.  (28),  the  pairings 


V  ’  iu0i 


»01 


~0i  *  ~0i 


i* 1)2 ^ 

V  ‘  -  ‘“oi  •  Soi'  "  ?oi  ■  ?oi'  *  ?oi  '  »oi  • 

i“l,2,. . . ,n 


(55) 


The  first-order  perturbation  solution  is  obtained  by  introducing 
-Oi  *  ~0i’  1  *  1*2, ...,2n  into  Eqs.  (41)  and  (43).  The  first-order 
eigenvalue  perturbation  is  . 


Xli  “  2  y0iAl~0i  *  1*1*2* • • • »2n  (56) 

H  _  t 

where  u^  =  u^^.  Because  is  real  and  symmetric,  all  A^  are  real. 
Furthermore,  assuming  that  the  matrix  C  is  positive  semidef inite,  A^  is 
negative  semldefinlte  and  hence  the  A^  are  nonpositive.  Moreover,  due 
to  the  real  symmetric  nature  of  the  matrix  A^,  we  can  restate  Eqs.  (56) 
as 


1  ha 

—  u„.A,u„. 

2  ~0i  l~0i 


Is 1,2,. . . , n 


which  is  to  say  that  the  first  order  perturbations  to  A^ 
are  the  same.  These  familiar  results  are  reassuring. 

The  coefficients  in  the  expansion  (38)  for  u^ 
u  HA  u 

eik  "  ?a~k  "K°S  »  i,k-l,2 . 2n,  i  i  k 


and  A 


are 


Oi’ 


(57) 


(58) 
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Recalling  the  symmetry  of  and  the  imaginary  nature  of  the  A^'s,  one 
can  show  that 


eik  "  eki  * 

i,k*l,2,.  • . ,2n 

(59a) 

ei'k'  "  eik 

t  ip k» 1 p 2  f •  •  0  ,n 

(59b) 

eik'  "  ei'k 

y  If  k* 1)2  y » • •  ,n 

(59c) 

Thus  if  we  regard  the  2n  x  2n  array  of  e's  as  partitioned  into  four  n  x  n 
arrays,  we  need  calculate  only  the  upper  triangular  elements  in  the  two 
upper  n  x  n  arrays.  Furthermore,  from  Eqs.  (58),  we  note  that  e  and 

IK 

eifc,  are  related  by  the  common  vector  A^u^.  Let  us  write  Eqs.  (38) 
as 

n 

~li  "  j^j^ilc-Ok  +  eik'-Ok^  *  *"l»2,...,n  (60a) 

Then,  we  can  show  that 

•  i“l»2,. . .  ,n  (60b) 

and,  from  Eqs.  (39),  (41)  and  (59a),  that 

yli  m  ,  1*1,2,. . . , 2n  (61) 

Next,  let  us  denote  the  unperturbed  eigenvectors  as  follows: 

~0i  “  ?0i  +  ±s0i  ’  ~0i'  “  ?0i  "  i?0i  ’  i"l»2,...,n  (62) 

T  H 

and  note  at  this  point  that  the  normalization  *  Yo^Ygi  "  ^  can 

realized  by  setting 

?0i?0i  "  *01*01  "  1  »  i-l,2,...,n  (63) 

Inserting  Eqs.  (62)  into  Eqs.  (43),  we  obtain 
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li.I 

-Oi 

?0i 

1  1 

f  I,  .  T. 

hv  2 

H  A1 
-Oi’ 

?0i’ 

2  I 

i?01Al?0i  +  201Al50iJ 

9 

,  i*1.2 . n  (64) 


Because  the  matrix  is  negative  semidefinite,  Eqs.  (64)  confirm  that 
all  first-order  perturbations  to  the  eigenvalues  are  real  and  nonpositive. 
Moreover,  inserting  Eqs.  (62)  into  Eqs.  (58),  we  can  write  for  1  i1  k 


-0kAl-0i  ~  ~0kAl-0i 


-  i 


~0kAl'0i  +  ~0kA3 


'ik 


2(“0i  ■  V 


1  <  i,  k  <  n  (65a) 


?okAi?oi  +  ?okAiy0i 


-  i 


'ik 


-0kAl~0i  ~  -0kAl-0i 


2(u0i  +  “W 


1  <  i  <  n,  n  +  1  <  k  <  2n  (65b) 


'ik 


f^0kAl-0i  +  -Ok^Oi  +  1  y0kAl?0i  ~  -0kAl-0i) 


2(<1>0i  +  V 


e 


-0kAl-0i  ~  -OkAl-Ql 


l<k<n,  n+l<i<2n  (65c) 


+  11^0kAl-0i  +  ~OkAl-Oi 


ik 


2(a,0i  *  “0k} 


n  +  1  <  i,  k  <  2n  (65d) 

Introducing  Eqs.  (62)  and  (65)  into  Eqs.  (38),  we  obtain 

rr 


Su"  z 

1  k-1 

k*i 


^0kAl-0i  "  -0kAl-0i  *  1f?0k*l?01  +  ~0kAl-0i  .  . 

- —  -  ■■  ' - ^  <?0k  +  ‘W 


2<“0i  *  "Ok’ 


+  (?OkAlS 


Z0i  *  -0kA1^0i 


-  i 


*0kAl*0i  ~  50kAl50i 


2(u  +  u  ) 

^  Oi  0k; 


■  (*0k  -  ‘W 

1*1,2,. . . ,n  (66a) 
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(66b) 


>li'  "  -li  *  i-l,2,...,n 

The  second-order  perturbation  eigensolutions  are  obtained  from  Eqs. 
(48),  (50)  and  (51).  From  Eqs.  (50),  (55)  and  (59a),  we  can  write 


X2i  '  “"“Ok  -  “01>eikeik 

K“1 


(u0k  +  “oi)eik*eikfl 


• » >n  (67a) 


X2i'  “  X2i  » 


(67b) 


It  should  be  evident  from  Eqs.  (62)  that  “  l,2,...,2n)  is  imaginary. 

From  Eqs.  (51)  and  (59)  we  can  write 


eii  "  Yii  “  2  ntf1[eimeim  +  eim,eim’]  * 


(68a) 


*ei'i'  *  Yi'i'  *  eii  *  i“l»2,...,n 


(68b) 


We  note  that  and  y±1  (i  ■  l,2,...,n)  are  real.  Comparing  Eqs.  (48a) 
and  (48b),  we  can  show  that 

®ik  *  Yik  ’  i.k-1,2,... ,2n,  i^k  (69] 

This  result,  coupled  with  Eqs.  (68)  ,  allows  us  to  state  that 


eik  "  Yik  *  1»kpl,2,.’.,2n 


and  hence  that 


^2^  "  ^21  *  i“l»2,« • • ,2n 


Let  us  now  write  Eq.  (48a)  as 


eik  "  *01-*0k  j(Xlk  "  Xll)6ik  +  J^<X0.  ~  X0k)eimemk 


"  <X0m  +  X0k)eim,em,k1 


i,k-l,2,. . .  ,2n,  i  k 


One  can  now  show  that 


“  eijc  *  i»kml,2,. . . ,n 


eik’  “  Vk  • 


(73a) 


(73b) 


From  Eqa.  (47a)  and  (70),  we  can  write,  in  analogy  with  Eqs.  (59a) 

—  n 

e,.  -  -  e,.  +  J  (e.  e  .  +  e.  ,e  ,.) 

ki  lk  ,  km  mi  km'  mi 

OF1 

i,k-l,2,...,2n,  i  A  k  (73c) 

Once  again,  if  we  regard  the  2n  *  2n  array  of  e's  as  partitioned  into 
four  n  x  n  arrays,  we  need  calculate  only  the  upper  triangular  elements 
in  the  upper  two  n  *  n  arrays.  Use  of  Eqs.  (73)  allow  us  to  fill  in  the 
remaining  elements.  Let  us  write  Eqs.  (45a)  as 
n 


~2i 


1  (*ik~0k  +  ^k'-Ok^  '  i-l,2,...,2n 


(74) 


k”l 


Bearing  in  mind  that  u^,  ■  u^,  Eqs.  (73a),  (73b)  and  (74)  allow  us  to 
state  that 


~2i'  "  ~2i  * 


(75) 


It  may  be  convenient  to  exhibit  the  real  and  imaginary  parts  of  the 
second-order  perturbations  explicitly.  To  this  end,  let  us  write 


eik  ”  ®Rik  +  leiik  »  . .  ,2n 


(76) 


where  subscripts  R  and  I  denote  the  respective  real  and  imaginary  parts 


of  eifc.  Substitution  of  Eqs.  (76)  into  Eqs.  (67)  yields 


X2i  "  1  ^  i“0k  "  V 


LeRik  +  eiikJ  + 


fok  ■  “oij 


2  2 
eRik’  +  *Iik 


•)] 


i*l,2 , . . . ,n  (77a) 
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(81b) 


“21*  *  *  i*l»2,...fn 

which  completes  the  second-order  perturbation  solution. 

It  should  be  observed  from  the  above  that  all  the  quantities  needed 
for  the  evaluation  of  the  first-  and  second-order  perturbations  are 
computed  by  means  of  the  eigenvalues  and  eigenvectors  of  the  unperturbed 
system,  i.e.,  those  of  the  undamped  gyroscopic  (or  nongyroscopic)  system. 


6.  System  Response  Based  on  the  Second-Order  Perturbation  Eigensolution 
In  Sec.  2,  we  obtained  the  equations 

t 


Zi(t)  -  eXitzt(0)  +  /0  eXi(t_T)Zu(T)dT  ,  i-1,2 . 


2n 


(82) 


It  will  prove  convenient  to  express  Eqs.  (82)  in  terms  of  the  state  vector 
and  the  associated  excitation  vector.  From  the  results  of  Sec.  2,  we  can 
write 

fq(o)' 


1  Tt  t 

zi(0)  "  2  YiL  Lq(o). 


1  T  -1 

Zu(t)  -  -  YiL 


Q(t) 


i-1,2,. ..,2n 

,  i-1,2, . . .  ,2n 

Equations  (82)  can  now  be  written  as 


z±(t) 


„  *it  1  T  T 

,^.r 

+  tC  eH(t-r)  I  ,T  -1 

Q(t) 

2  -i 

Lg(o)J 

j0  2  -i 

dr 

i-1,2,... ,2n 


(83) 

(84) 


(85) 


which  permits  us  to  write  the  state  vector 
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|Ct> 

Ls(t) 


,-T  1  r-T  2v  i  Aic  TtT 

L  E  u.z  (t)  «  ■=•  L  E  ie  u.v  L 

i-1  1  z  i-1 


q(0) 

,q(o) 


Ai(t-t)  T-l 
e  i  u^v^L 


+  ^0 

From  Eq .  ( 36a) ,  we  have 


Q(t) 

0 


dx 


(86) 


Xi  "  X0i  +  Xli  +  X2i  +  •**  *  i”1,2’* ’ ’ *2n 


(87) 


In  Sec.  5,  we  noted  that  A  and  A^  are  imaginary  quantities,  while  A^ 
is  real  and  negative.  Furthermore,  we  noted  that 


0i’ 

“  X0i  *  ‘  X0i 

,  i«l,2,...,n 

li’ 

*  Xli 

,  i- 1,2,. . . ,n 

(88) 

2i’ 

"  X2i  "  “  X2i 

,  i-1, 2 , • • • ,n 

where  the  meaning  of  a  primed  subscript  was  indicated  in  the  previous 
section.  In  view  of  this,  let  us  express  Eqs.  (87)  as 


Ai  *  “  Yi  "  ia,di  ’ 

V  »  \  ”  -  y±  +  i“di  .  i=-l,2,...,n 


(89) 


where 

y±  ~  ~  Xli  >  i-l,2,...,n 

(90) 

iudi  "  X0i  +  X2i  *  i"1.2.-.-.n 

and  we  note  that  the  a),,  represent  frequencies  of  damped  oscillation. 

al 

From  Eqs.  (36b),  we  have 

“i  *  y?0i  +  -li  +  -2i  *  i“l»2».  •  •  »2n  (91) 

where  we  have  shown  that 


2k 


-Oi'  "  ?0i 

,  i* 1,2,...  ,n 

-li'  "  -li 

,  1*1*2 , ...  ,n 

(92) 

~2i'  “  ~2i 

,  i* 1,2,...  ,n 

we  can  state 

that 

Si»  =  “i  * 

i*l ,2,... , n 

(93) 

Furthermore,  let  us  write 
2n 

-i  *  J  Eik-0k  ’  i=l»2,...,2n 

k=l 


where 


Eik  *  6ik  +  eik  +  eik  +  ’  i*W,2»"’»2n 


From  Eqs.  (36c) ,  or 


‘  ?01  +  ?U  +  ?2i  +  •••  ■  1*1’2 . 2" 


we  can  similarly  show  that 

Y1*  *  Yi  ,  i“l»2, . . . ,n 

and  we  can  write 


?1  ■  2,°ltSok  -  1-1'2 . 2” 

k“l 


where 

s  6^  —  6^  +  +  ...  ,  i,k“l,2,. . .  ,2n  (99 

Let  us  now  reconsider  Eq.  (86).  From  Eqs.  (89),  (93)  and  (97)  we 
conclude  that  the  i'th  term  is  the  complex  conjugate  of  the  ith  term, 
for  every  i  and  i'.  Hence,  we  can  rewrite  Eq.  (86)  as 
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(t)  T-TR  "  j  Xit  I  I  ?(0) 
-  =  L  Re  S  ■je  1  u  v  L  - 

(t)J  i-ll  _q(0) 


r  X^t-t)  T  -1  ?(T) 

+  I0  e  Vi1  —  < 

L  o 


From  Eqs.  (89),  we  can  write 

e*iC  =  e  Yit(cos  o>dit  -  i  sin  u>dit)  ,  i*l,2,...,n 


and  from  Eqs.  (94)  and  (98),  we  have 


(100) 


(101) 


I  E..  urtl 
k«l  ik'0k 


2n  —  —  T 

2  Gim-0m 
m=l 


1-1,2,. 


(102) 


so  that,  using  Eqs.  (62),  we  can  express  Eqs.  (102)  as 


T  n  r  f  V 

'iYi  "  k,^iLEikl-°k  +  i?okJ  +  Eikl-°k  '  i?0kL 

■  -  *ioi)  +  (?oI  +  i5oi]]  - 


i-1,2,. . . ,2n 


(103) 


Bearing  in  mind  that  the  coefficients  E^  and  G  are  complex  quantities. 


we  can  write 


uy^  «  +  il^  ,  i»l,2,...,2n 


(104) 


where  the  2n  *  2n  matrices  R^  and  1^  are  the  respective  real  and 

T 

imaginary  parts  of  the  2n  x  2n  matrix  u^.  Introducing  Eqs.  (101)  and 
(104)  into  Eq.  (100),  we  obtain 
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q(t) 

,q(t)J 


_  n 
— T 

L  Z 
i-1 


e  Y*t  cos  w .  t  R.LT 
Ql  1 


+  e  Yit  sin  oj,.t  I.L 

Ql  1 


+  /0  e  Yi^c  T^cos  a)dl(t-T)RiL_1 
+  /Q  e"Yl(t"T)sin  “^(t-x)!^"1 


'q(O) 
S(°)J 
q(o) 

LsC0)j 

Q(t) 


o 


dx 


dx 


(105) 


which  expresses  the  response  in  terms  or  real  quantities  alone. 

Note  that  the  above  formulation  remains  valid  in  the  case  of  non- 
gyroscopic  systems,  G  *  0. 


7.  Numerical  Examples 

Consider  the  two-degree  of  freedom,  damped  gyroscopic  system  of 
Ref*  1  (shown  here  in  Fig.  1).  The  equations  of  motion  can  be  written 
in  matrix  form  as 


Mq(t)  +  (G  +  C)q(t)  +  Kq(t)  *  Q(t) 


(106) 


where 


M 


m  0 
0  m 


C 


c  0 

0  0. 


G 


K 


’  0  -2mri 

_2m0  0 

kj-mO2  0 
0 


(107) 


Let  us  take  m  -  1kg,  ft  -  lrad  s-1,  c  *  0.1kg  s-1,  ^  *  3kg  s“2,  k2  ** 
_2 

4kg  s  ,  and  note  that  k^,k2  and  0  have  been  picked  such  that  matrix 
K  is  positive  definite.  From  Eqs.  (3)  we  form  the  matrices  M*  and  K* 
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Computations  of  the  perturbations  to  the  undamped  elgensolutlons 


proceeds  as  Indicated  previously.  The  results  are  summarized  In  Table  1. 
One  should  note  that  the  system  eigenvalues  and  associated  eigenvectors 
occur  In  complex  conjugate  pairs.  Hence,  we  Include  only  half  of  them. 
Furthermore,  we  Include  only  the  right  eigenvectors. 

Inspection  of  Table  1  reveals  the  quality  of  the  convergence.  It 
Is  also  Interesting  to  Inspect  convergence  via  a  biorthonormality  check 

of  the  perturbed  eigenvectors.  In  Table  2,  we  present  the  matrix 

1  T 

2  V  U  .  Numerical  values  within  parentheses  represent  products 

2  (y0j  +  ylj  +  y2j) (»0i  +  yli  +  y2i)  '  . 2n  (llla) 


while  values  not  within  parentheses  represent  products 
2  (Y0j  +  ylj)  (v0i  +  -li)  »  . 2n 


(111b) 


The  response  to  excitation  in  the  form  of  the  Dirac  delta  function, 
with  Impulse  equal  to  1  kgms  2  has  been  computed.  Figure  2  represents  a 
plot  of  the  coordinate  y(t)  versus  t  due  to  an  excitation  applied  In  the 
y-direction.  The  two  curves  represent  the  response  obtained  by  the  gen¬ 
eral  theory  and  the  0(0)  response.  Within  the  accuracy  of  this  plot,  the 
0(0)  +  0(1)  and  0(0)  +  0(1)  +  0(2)  'responses  are  Identical  to  that  obtained 
by  the  general  theory. 

As  another  example,  let  us  consider  Che  slightly  damped  nongyroscopic 
system  depicted  in  Fig.  3.  The  equations  of  motion  considered  are 


Mq(t)  +  Cq(t)  +  Kq(t)  -  Q(t) 


(112) 


where 
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i  o]  r  0.4  -0.2]  r  5  -4 

M  -  ,  C  -  ,  K  -  (113) 

0  2j  L-0.2  0.2J  j_-4  4 


Matrices  Aq  and  are 


Summaries  of  the  exact  and  perturbed  elgensolutions  are  given  in  Table  3 
and  a  biorthonormality  check  is  given  in  Table  4.  The  quantities  displayed 
in  this  table  have  the  same  meaning  as  those  in  Table  2.  The  response  of  the 
system  to  an  excitation  in  the  form  F^(t)  *  0,  F^(t)  *  6(t)  is  presented 
in  Fig.  4.  As  in  Fig.  2,  the  response  obtained  by  the  general  theory  and 
the  0(0)  +  0(1)  and  0(0)  +  0(1)  +  0(2)  responses  are  all  given  by  one 
curve,  whereas  the  0(0)  response  by  the  other  curve. 


Table 


15.521 


Table 


o 

NO 

*4 

co9  1 

CM 

© 

CM 

oo 

O' 

NO 

O' 

m 

ON 

O' 

o 

CO 

CO 

so 

CO 

a'>aoto-a,cNi/ir~cN 

>»0<NI,1NO'» 

oonor'H^av 

A  •  e  •  e  i  *  • 

vO©r*00\©fH© 

m  l  oo  i  ©  +  ©  + 


'J  O  H  o  H  n  n 

O'  CM  ©  *4  O'  CM 

o  o  n  o  O'  O'  o 

*  •  •  •  •  I  • 

COOmOHOH 
ifl  *r4  H  ■H  "H  H 

oo  i  »n  10  +  0 


NO 

CM 

00'  ' 

CM 

© 

oo 

ON 

*4 

to 

in 

ON 

co 

co 

o 

o 

NO 

rv. 

© 

tj\  00  0V®O«a-CMlAf>«c4  o  O  O' 

O  00  OiOO'NHNO'O'  u'  ff. 

vO  u">  <IOOOONrlsro\  OO  O  4 

M  »  #»•*•»•*•  At* 

cm  cm  mO'COO'OHo  so  ©  m 

CM  «H  H*Hin«HCM«HCM*H  lA  *H  H 

cm  i  m  i  oo  io  +  o+  oo  i  m 


On  ©  CM 

on  ©  sr 

nr  co  on 

«  e  * 

m  no  m 

rH  m  CM 

m  oo  »H 


©  on  00  CM 

on  O'  *  n 

oo  <r  o  on 

d  *  »  * 

no  m  cm  m 

m  rH  ON  CM 

oo  m  on  rH 


m 

n* 

m 

ON 

m 

00 

CM 

© 

«» 

«* 

* 

•4  m 

NO  00 

rH  00 

©  -4 

©  r* 

m  © 

■4  m 

m  m 

©  © 

•  • 

*  • 

*  • 

n*  o 

CM  CM 

NO  © 

CM  *H 

CM  *H 

rH  *H 

CM  On  CM  CM  ON 

fo  o  n  h  m  o' 


I  ©  +  ©  + 


'  ©  00  O' 

O  00  ON 

oo  r-*  co 

«  *  • 

cm  O  ©  ©  rv  cn 

©  m  cm  ©  oo  O' 

h  ©  co  o  ©  on 

•  •  •  •  •  t 

oo  ©  m  ©  cm  © 

m  *H  H  yH  ,*H 

oo  i  m  i  o  V  P 


33 


-10.125,4581  I  II  -iO. 125,9421  -10.125,459 


(A.q  +  A^)u^  "  *  i“l»2,...,2n 

where  matrices  Aq  and  A^  are  real  and  where  A^  is  one  order  of  magnitude 
smaller  than  A^.  The  analysis  was  based  upon  knowledge  of  the  eigen- 
solutions  when  is  the  null  matrix.  The  perturbation  theory  was  applied 
to  slightly  damped  gyroscopic  systems.  In  this  case,  A^  is  skew 
symmetric  and  A^  is  symmetric,  so  that  special  computational  advantages 
can  be  realized.  Note  also  that  the  nongyroscopic  systems  can  be  handled 
within  the  context  of  the  same  general  formulation  by  simply  letting 
G  =  0. 

As  an  example,  a  two-degree  of  freedom  slightly  damped  gyroscopic 
system  was  analyzed.  As  another  example,  a  nongyroscopic  system  has  been 
treated.  Even  for  relatively  large  damping,  the  perturbation  results  for 
both  cases  agree  well  with  the  solutions  obtained  by  algorithms  for 
general  matrices.  Because  the  present  formulation  is  based  on  the  eigen- 
solution  for  real  symmetric  matrices,  it  should  prove  far  superior  for 
high-order  systems. 
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Figure  1.  Gyroscopic  System 

Figure  2.  Response  of  the  Gyroscopic  System 

Figure  3.  Nongyroscopic  System 

Figure  4.  Response  of  the  Nongyroscopic  System 


FIGURE  4 


